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ABSTRACT 

A  significant  improvement  has  been  developed  in  the  first  few  months  of 
this  program  that  allows  the  phase  closure  technique  to  deal  more  correctly  with 
the  case  of  distributed  inhomogeneities.  Our  modification  of  the  phase  closure 
technique,  to  which  we  refer  as  "cross-correlated  subarrays"  involves  processing 
the  hydrophone  signals  in  groups,  corresponding  to  operating  the  array  as  a 
series  of  smaller  subarrays.  The  beamwidth  of  a  subarray  is  chosen  to  restrict 
the  response  of  the  system  to  the  angular  field  over  which  the  assumptions  of 
phase  closure  are  valid.  If  we  let  be  the  coherence  scale  of  the  system  (a 
function  of  oceanic  parameters,  together  with  range  R  and  wavelength  A.),  then 
for  phase  closure  to  be  valid  (with  or  without  cross-correlated  subarrays)  we 
require  Ax^  <  (XR)^/^.  Cross-correlated  subarrays  are  required  explicitly  if 
the  maximum  source  separation,  s,  is  such  that  s  >  ax^..  These  expectations  have 
been  tested  by  means  of  a  series  of  numerical  experiments.  The  results  are 
consistent  with  the  predictions,  giving  us  confidence  that  the  main  features  of 
our  algorithm  are  correct.  The  next  step  is  to  apply  the  algorithm  to  actual 
sonar  patrol  data,  to  be  obtained  from  Autonetics  Marine  Systems  Division  of 
Rockwell  International.  Pending  the  outcome  of  these  tests,  we  expect  that  the 
algorithm  should  provide  a  practical  and  effective  method  for  improving  the 
quality  of  passive  sonar  images.  , 
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INTRODUCTION 

The  purpose  of  this  program  is  to  explore  the  possibility  of  applying 
the  phase  closure  principle  to  passive  sonar  imaging  in  order  to  overcome  the 
distorting  effects  caused  by  both  the  propagation  medium  and  phase  errors  in  the 
instrument  itself.  The  phase  closure  technique  was  originally  developed  for 
Very  Long  Baseline  radio  astronomy  in  order  to  overcome  the  phase  distortion 
caused  by  the  Earth's  atmosphere  (Readhead  and  Wilkinson,  1978).  An  important 
difference  between  this  case  and  passive  sonar  is  that  for  radio  astronomy,  the 
distorting  medium  represents  a  thin  sheet  at  one  end  of  the  propagation  path, 
whereas  for  the  sonar  case,  the  inhomogeneities  responsible  for  phase  distortion 
are  distributed  continuously  between  the  source  and  receiver. 

Prior  to  the  beginning  of  this  program,  Rockwell  International  Science 
Center  had,  under  the  auspices  of  its  Independent  Research  and  Development 
funds,  already  laid  some  of  the  foundation  for  the  application  of  the  phase 
closure  technique  to  passive  acoustic  imaging.  In  this  work  (Marsh,  Richardson, 
and  Martin  1985)  the  theoretical  limits  of  the  standard  phase  closure  technique 
were  determined,  and  these  limits  demonstrated  by  performing  a  series  of  numeri¬ 
cal  experiments.  A  closely  related  technique  for  overcoming  phase  distortion 
effects  has  recently  been  discussed  by  Paul  raj  and  Kailath  (1985).  Their  tech¬ 
nique  models  the  phase  errors  in  exactly  the  same  way  as  the  phase  closure  tech¬ 
nique,  although  the  solution  to  the  imaging  problem  is  accomplished  by  eigenvec¬ 
tor  methods.  Also  closely  related  is  the  concept  of  phase  recovery  from  bi spec¬ 
tra  (Bartelt,  Lohmann,  and  Wirnitzer,  1984).  As  in  the  case  of  the  conventional 
phase  closure  method,  both  of  these  methods  would  be  applicable  to  instrumental 
phase  errors,  but  would  not  be  able  to  handle  the  case  of  distributed  inhomo¬ 
geneities  without  some  modification. 

Our  efforts  since  the  beginning  of  the  present  DARPA-funded  project 
have  been  concerned  with  extending  the  validity  of  the  phase  closure  technique 
to  the  case  of  distributed  inhomogeneities.  Our  principal  accomplishments  have 
been  : 
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1.  Extend  the  range  of  validity  of  the  phase  closure  technique  using 
cross-correlated  subarrays. 

2.  Theoretical  understanding  of  the  behavior  of  the  algorithm  in  the 
context  of  distributed  inhomogeneities. 

3.  Numerical  simulation  of  the  "cross-correlated  subarrays" 
al  gorithm. 

We  will  discuss  the  theoretical  basis  of  this  work,  and  present  the  results  of 
numerical  testing  with  synthetic  data. 


THEORETICAL  CONSIDERATIONS 


As  discussed  by  Marsh,  Richardson  and  Martin  (1985),  the  far-field 
intensity  distribution  I(^)  of  a  spatially  incoherent  sound  source  in  the  case 
of  an  inhomogeneous  ocean  is  related  to  the  measured  visibility  by: 


^meas 

ij 


*  -2niu--.e„ 

~  I  9-in  ® 

^  ^in  ^jn  ~n  n 
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where  e^  is  a  unit  vector  in  the  direction  of  the  nth  source  component  as  viewed 
from  the  receiver,  is  the  solid  angle  subtended  by  the  nth  source,  u^-  is 
the  baseline  vector  joining  hydrophones  i  and  j  in  units  of  wavelengths,  and  g^p 
and  gjp  are  complex  gain  factors  representing  the  effect  of  propagation 
anomalies  in  the  medium. 


The  principal  assumption  underlying  the  standard  phase  closure  tech¬ 


nique  is: 


^in  “  ^i 
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i.e.,  the  gain  factors  vary  more  slowly  with  source  position  than  with  receiver 
position.  If  we  assume  that  the  principal  effect  of  the  medium  is  a  phase 
distortion,  i .e., 

id.. 

Qi  =  e  '  (3) 

then  Eq.  (2)  implies  that  we  can  express  the  measured  phase  =  arg(v'l’^^^) 

tPUG 

in  terms  of  the  true  phase  (corresponding  to  a  homogeneous  ocean)  as: 

♦Tf'  =  4?' "  -  *0  •  <«) 

The  standard  phase  closure  technique  attempts  to  obtain  the  phase  errors  4.^ 
using  a  measurement  model  of  the  form  (4),  and  correct  the  data  accordingly. 

For  a  given  set  of  physical  conditions,  one  can  define  an  angular 
field-of-view  over  which  (2)  is  satisfied.  Provided  all  of  the  acoustic 
sources  are  located  within  this  narrow  field,  the  standard  phase  closure  tech¬ 
nique  is  valid.  Unfortunately,  this  is  unlikely  to  be  true  in  practice,  since 
sources  occur  typically  over  a  wide  range  of  azimuths.  Recent  work  has  been 
devoted  to  the  problem  of  extending  the  phase  closure  technique  to  handle  this 
situation.  A  conceptually  simple  way  of  accomplishing  this  is  to  divide  the 
array  up  into  a  number  of  identical  subarrays,  and  operate  each  as  a  small 
phased  array,  with  the  beams  all  pointed  in  the  same  direction.  The  beam 
pattern  of  each  subarray  will,  of  course,  be  much  broader  than  the  resolution  of 
the  full  array.  We  will  refer  to  the  beam  pattern  of  each  subarray  as  the 
primary  beam.  We  will  treat  each  subarray  as  an  individual  element  in  a  larger 
scale  array,  and  cross-correlate  the  signals  from  each  subarray  as  if  we  were 
dealing  with  an  array  of  single  hydrophones.  The  effect  of  this  is  that  the 
angular  response  of  the  larger  array  will  have  been  multiplied  by  the  primary 
beam  pattern  corresponding  to  an  individual  subarray.  By  an  appropriate  choice 
of  the  number  of  elements  in  a  subarray,  we  could,  provided  certain  conditions 
are  met,  restrict  the  primary  beamwidth  A0pp^-  to  the  field  AGp^.  over  which  (2) 
is  satisfied.  In  the  case  of  a  homogeneous  ocean,  A0pp.j  can  be  obtained  simply 
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from  the  Fourier  transform  of  the  aperture  distribution  for  a  single  subarray. 
In  the  case  of  an  inhomogeneous  ocean,  however,  the  phase  distortions  will 
broaden  the  response,  ultimately  setting  a  lower  limit  A9(j  on  the  primary  beam- 
width  obtainable. 

Both  A9pj.  and  A9|j  can  be  related  to  the  phase  structure  function 
D(x  -  x'),  defined  by  Flatte  et  al .  (1979)  as: 

D(x  -  x')  =  E[(i.(x)  -  4;(x')]^  (5) 


where  x  and  x'  represent  the  positions  of  two  receivers  located  at  a  distance  R 
from  a  point  source  and  displaced  perpendicular  to  the  direction  of  propagation, 
(i-(x)  represents  the  spatial  variation  of  the  phase  of  the  received  signal  and  E 
is  the  expectation  operator.  We  can  then  express  A9pj.  and  A9|3  in  terms  of  some 
convenient  criterion  limiting  the  amount  of  permissible  phase  variation  Ad;  over 
the  source  and  receiver,  respectively.  Provided  the  inhomogeneities  have  the 
same  statistical  character  throughout  the  region  between  source  and  receiver, 
and  taking  Ad;  to  be  1  radian,  we  can  then  write 


R  A9 


pc 


\ 

A9b 


AX 


C 


where  \  is  the  wavelength  and  ax^.  is  defined  by 


(6) 


D(AxJ  =  1  .  (7) 

Provided  A9pj.  >  A9^,  the  primary  beamwidth  can  be  restricted  to  A9p^,  enabling 

the  phase  closure  technique  to  produce  an  image  over  the  restricted  field 

9p  -  A9pj./2  to  9p  +  A9pj./2  where  9p  represents  the  center  of  the  restricted 

field.  A  mosaic  of  subimages  could  then  be  produced,  corresponding  to  various 

values  of  9p,  and  the  final  image  produced  by  a  combination  of  these  subimages. 

In  order  to  produce  the  r^^  subimage,  the  phase  closure  technique  would  be 

( rl 

applied  to  the  set  of  visibilities  V:.'  defined 

*  w 


5 

C7173A/jbs 


Rockwell  International 

Science  Center 


SC  5427. IR 


=  I  I  a^(i)  a  (j)  E(S.S  ")  e 

^  k=k  (i)+l  Jl=JL(j)+l  ^  ^  ^ 


k„{i)+ns 


-2«UkjSln  6^ 


where 

ko(i)  =  (i  -  Dn^ 

Ao(j)  =  (j  -  Dn^ 

and  S|^,  represent  the  signals  received  by  the  and  hydrophones;  a|(,(i), 
ajj^(j)  are  weights  representing  the  apodizing  function  for  each  subarray,  n^  is 
the  number  of  elements  in  each  subarray. 

The  new  algorithm  sketched  above  would  enable  phase-closure  imaging 
over  a  wide  field-of-view,  and  should  provide  a  substantial  improvement  in 
imaging  performance  over  that  obtained  by  conventional  beamforming.  We  will 
refer  to  it  subsequently  as  "phase  closure  with  cross-correlated  subarrays."  We 
now  consider  the  question  as  to  the  physical  regime  over  which  this  algorithm 
will  be  valid.  A  useful  quantity  to  bear  in  mind  in  this  regard  is  defined 
by  (7),  which  represents  the  coherence  scale  for  the  system  at  the  particular 
range  and  frequency.  The  condition  that  A0p(-  >  Ae^,  is  equivalent  to 

Ax^  >  .  (9) 

Equation  (9)  thus  represents  a  criterion  by  which  one  can  determine 
whether  or  not  phase  closure  (with  cross-correlated  subarrays)  is  capable  of 
improving  a  distorted  image.  In  order  to  relate  it  to  more  fundamental  physical 
parameters,  however,  requires  a  knowledge  of  the  form  of  the  phase  structure 
function  D,  which  in  turn  is  related  to  the  spatial  correlation  function  of  the 
refractive  index.  If  we  assume  that 

E(i,(x)  (i,(x')  =  exp(-(x-x' )^/L^)  (10) 
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(11) 


and  hence 

AX  =— L  .  (12) 

If  the  inhomogeneities  are  statistically  uniform  and  isotropic,  then 
from  Flatte  et  al.,  (1979): 

$2  ,  0.4  <u^>  RL  (13) 

where  >  represents  the  variance  of  the  refractive  index  fluctuations. 
Equations  (12)  and  (13)  then  give 

AX  =  0.2  \  (— b - )^^^  .  (14) 

<  />  R 

Thus  from  (9),  we  find  that  the  maximum  range  over  which  the  phase 
closure  technique  (with  cross-correlated  subarrays)  can  be  applied  to  a  distri¬ 
buted  medium  is 


Rmav  =  0-2 
max 


(•^) 

> 


1/2 


(15) 


In  the  case  of  horizontal  imaging  in  the  ocean  at  500  Hz  (assuming 
L  =  5  km  and  <jx2>l/2  =  5  x  10“^)  we  obtain  R^gj^  =  49  km,  which  suggests  that 
substantial  improvement  in  sonar  images  from  horizontal  arrays  should  be 
possible  over  ranges  up  to  this  value. 


& 


V'-.  ’ 


AV- 
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NUMERICAL  EXPERIMENTS 

In  order  to  test  the  concepts  discussed  in  the  previous  section,  a 
series  of  numerical  experiments  was  performed.  In  these  experiments,  an  assumed 
source  consisting  of  4  pointlike  components  spread  over  a  distance  of  200  m  was 
observed  at  a  distance  of  1  km  and  a  frequency  of  300  Hz,  using  a  regularly 
spaced  array  of  length  200  m  containing  21  hydrophones.  The  sound  velocity  was 
assumed  to  be  1500  m/s.  Ocean  inhomogeneities  were  simulated  by  placing  four 
equally  spaced  phase  screens  between  the  source  and  receivers.  Along  each  phase 
screen,  the  phase  was  made  to  vary  randomly,  with  a  correlation  length  L.  The 
rms  phase  deviation  through  the  entire  system  was  chosen  to  be  approximately  1 
radian,  so  that  L  corresponds  to  Ax^.  Various  values  of  L  were  assumed,  in  some 
cases  chosen  to  purposely  violate  the  assumptions  for  phase  closure.  The 
propagation  was  calculated  according  to  the  Fresnel  approximation,  as  discussed 
by  Marsh,  Richardson  and  Martin  (1985).  In  the  case  of  the  cross-correlated 
subarrays  algorithm,  each  subarray  had  5  elements,  weighted  with  a  rectangular 
function,  i.e.,  no  apodizing  was  used,  A  mosaic  of  3  subimages  was  produced  in 
each  case,  spaced  by  the  half-power  widths  of  the  corresponding  primary  beams 
(sine  functions).  Simple  addition  was  used  in  order  to  construct  the  final 
image  from  the  3  subimages. 

Since  in  these  experiments  L  =  ax^,  Eq.  (9)  implies  that  phase  closure 
would  be  invalid  if  L  <  (XR)^/^.  In  the  case  of  L  >  (\R)^/^,  phase  closure  would 
be  applicable,  but  cross-correlated  subarrays  would  be  required  if  the  source  com¬ 
ponents  are  separated  by  a  distance  s  greater  than  L.  These  two  criteria  form  a 
rather  crucial  test  of  the  theoretical  basis  of  our  algorithm,  and  values  of  L 
were  chosen  to  test  them.  In  these  experiments,  (XR)^/^  corresponds  to  71  m. 

We  now  discuss  the  results. 

Case  a;  L  >  (XR)^/^,  s  =  L. 

The  assumed  value  of  L  was  200  m.  In  this  case  we  expect  phase  closure 
to  be  valid,  and  cross-correlated  subarrays  should  not  be  necessary.  The 
results  are  shown  in  Fig.  1,  which  fulfills  our  expectations.  The  use  of  phase 
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Case  b:  L  >  s  =  2L. 

The  assumed  value  of  L  was  100  m.  We  expect  phase  closure  to  be  valid, 
but  the  sources  are  now  so  widely  spaced  in  comparison  to  the  length  scale  of 
the  inhomogeneities  that  cross-correlated  subarrays  should  be  required.  The 
results,  shown  in  Fig.  2  support  these  expectations.  Standard  phase  closure 
gave  a  very  poor  result,  whereas  cross-correlated  subarrays  gave  a  reasonable 


reconstruction. 


SC86-32702 


TRANSVERSE  DISTANCE  (m) 

Fig.  2  Imaging  results  for  the  case  of  moderate  phase  distortion.  Parameters 
are  the  same  as  for  Fig.  1,  except  that  s  =  2L,  i.e.,  the  sources  are 
now  so  widely  spaced  in  comparison  to  the  length  scale  of  the  inhomo¬ 
geneities  that  cross-correlated  subarrays  are  required. 
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Case  c:  L  < 

The  assumed  value  of  L  was  50  m.  In  this  case  we  expect  phase  closure 
to  be  completely  invalid,  and  the  results  in  Fig.  3  show  this  to  be  true.  It  is 
interesting,  however,  that  the  cross-correlated  subarrays  image  bears  some 
resemblance  to  the  assumed  source,  although  the  separation  between  components  is 
incorrect. 

SC8S-32704 
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DISCUSSION  AND  FUTURE  PLANS 

The  numerical  simulations  discussed  above  give  us  confidence  that  the 
main  features  of  our  algorithm  are  correct,  and  that  substantial  image  improve¬ 
ment  should  be  possible  in  the  case  of  towed  arrays.  The  next  step  in  the  pro¬ 
gram  is  to  apply  the  "phase  closure  +  cross-correlated  subarrays"  algorithm  to 
real  data,  taken  from  sonar  patrols.  Data  for  this  purpose  exists  at  Autonetics 
Marine  Systems  Division  (AMSD)  of  Rockwell  International.  Initial  tests  will  be 
performed  on  the  assumption  that  sources  are  at  a  distance  of  the  order  of  10  km. 
For  typical  ocean  parameters  obtained  from  Flatte  et  al.,  (1979)  =  5  x  10"^, 

L  =  5  km) ,  at  a  frequency  of  500  Hz  we  obtain  ax^.  =  760  m,  corresponding  to  a 
subtended  angle  at  the  receivers  of  4.4°.  This  means  that  hydrophone  signals 
should  be  averaged  so  as  to  give  a  primary  beam  of  this  width.  If  the  hydro¬ 
phones  are  spaced  at  half-wavelength  intervals  (1.5  m),  this  means  processing 
the  hydrophone  signals  in  groups  of  about  20. 

Physical  parameters  corresponding  to  patrol  data  in  hand  will  be 
obtained  from  AMSD  in  the  near  future.  A  numerical  simulation  using  these 
parameters  will  be  performed  prior  to  testing  with  the  patrol  data. 
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